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Abstract 

We present a theory describing the superconducting (SC) interaction of Dirac elec- 
trons in a quasi-two-dimensional system consisting of a stack of N planes. The 
occurrence of a SC phase is investigated both at T = and T 7^ 0, in the case of 
a local interaction, when the theory must be renormalized and also in the situation 
where a natural physical cutoff is present in the system. In both cases, at T = 0, we 
find a quantum phase transition connecting the normal and SC phases at a certain 
critical coupling. The phase structure is shown to be robust against quantum fluc- 
tuations. The SC gap is determined for T = and T 7^ 0, both with and without 
a physical cutoff and the interplay between the gap and the SC order parameter 
is discussed. Our theory qualitatively reproduces the SC phase transition occurring 
in the underdoped regime of the high-Tc cuprates. This fact points to the possible 
relevance of Dirac electrons in the mechanism of high-Tc superconductivity. 
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1 Introduction 



Surprisingly, there are many condensed matter systems in one and two spatial 
dimensions containing electrons that may be described by a relativistic, Dirac- 
type lagrangian, namely Dirac electrons. Even though these are evidently non- 
relativistic systems, this fact occurs because in some materials there are special 
points in the Brillouin zone where two bands touch in a single point around 
which the electron dispersion relation may be linearized as e{k) = vp\k\. The 
kinematics of these electrons can be described by a Dirac-type lagrangian 
where the velocity vp determines the angle of the Dirac cone. At the tip of 
this cone the Fermi surface reduces to a point, the Fermi point, and the density 
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of states vanishes. The elementary excitations around a Dirac point are Dirac 
electrons. They are, after all, a result of the electron-lattice interaction. 

There are many important quasi-two-dimensional systems containing Dirac 
electrons. Among them wc could mention the high-Tc cuprates, where Dirac 
points appear in the intersection of the nodes of the d-wave superconduct- 
ing gap and the 2D-Fcrmi surface. Because of these nodes the low-energy 
quasiparticle spectrum is gapless and the dispersion relation can be linearized 
as described above [1]. The low-energy quasiparticle dynamics is determined 
exclusively by these points, since they are occupied even at very low tem- 
peratures [2,3,4,5]. For these reasons, Dirac electrons are expected to play an 
important role in the cuprates. This motivates the application of the model 
introduced below to these materials and the results thereby obtained seem to 
point towards this direction. 

Dirac electrons also appear in semi-metals such as graphene sheets or stacks 
thereof, namely graphite, where the vanishing density of states at the Fermi 
points has important consequences in the electronic properties like, for in- 
stance, the absence of screening of the Coulomb potential [6,7,8,9,10,11]. 

Another class of materials where the presence of Dirac electrons produces 
interesting effects are the rare-earth dichalcogenides such as 2H—TaSe2, 2H— 
NhSe2, 2H — TaS2 and 2H — NbS2. In these systems charge-density- wave 
order coexists with superconductivity at low temperatures and Dirac points 
form in the intersection of the Fermi surface with the nodal lines of the charge- 
density- wave order parameter [12,13]. This is a particular example of nodal 
liquids, which in general contain Dirac electrons [14]. Finally, we may list 
carbon nanotubes as another important type of materials that have also been 
shown to possess Dirac electrons [15]. 

In this work, we present a theory describing the superconducting interaction 
of Dirac electrons associated to two distinct Dirac points belonging to a stack 
of N planes. We analyze the conditions for the existence of a superconducting 
gap, both at T = and T ^ 0, either for finite N or in the limit N ^ oo, 
which corresponds to the case when the system is actually three-dimensional. 

At T = 0, we show that the system presents a quantum critical point sepa- 
rating the normal and superconducting phases and determine the supercon- 
ducting g function of the coupling constant. The theory is renormal- 
ized in a 1/A^ expansion. A renormalization group analysis is then performed, 
demonstrating the independence of physical quantities from the renormaliza- 
tion point. We also investigate the effect of quantum fluctuations in our results 
and demonstrate that the phase diagram found in mean field is robust against 
these fluctuations. 

We then consider the finite temperature case and determine the superconduct- 
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ing gap A as a function of the temperature and of the zero temperature gap. 
We find a critical temperature Tc, where the gap vanishes. The possibihty of 
occurrence of dynamical generation of a superconducting gap without the cor- 
responding spontaneous breakdown of the U(l) symmetry, in compliance with 
the Coleman-Mermin-Wagner-Hohenberg [16] theorem is discussed in detail, 
as well as the associated Kosterlitz-Thouless transition suffered by the phase 
of the complex order parameter. 

We finally make a detailed study of the situation usually found in condensed 
matter applications, where a natural momentum cutoff exists in the system, 
both at T = and T 7^ 0. We consider the general regime, where the cutoff is 
not necessarily much larger than the gap, as well as the weak coupling regime 
where the cutoff is much larger than the gap. The former situation is likely to 
be relevant for high-Tc superconductors while the latter is the one found in 
conventional BCS superconductors. 

The quantum phase transition occurring in our model and the behavior of Tc 
around the quantum critical point qualitatively reproduce very well the super- 
conducting transition in the high-Tc cuprates in the underdoped region. This 
suggests that Dirac electrons may play an important role in the mechanism of 
high-Tc superconductivity. 



2 The Model and the Effective Action 



We consider a quasi- two- dimensional electronic system consisting of a stack 
of planes containing two Dirac points. In addition, we introduce an internal 
index a = 1, A^, supposed to characterize the different planes to which the 
electrons may belong. The electron creation operator, therefore, is given by 
i/jj„a, where i — 1,2 are the Dirac indices, corresponding to the two Fermi 
points, o" =1, 1, specifies the z-component of the electron spin and a = 1, N 
labels the electron plane. We assume, further, that there is a BCS-type super- 
conducting interaction, whose origin is understood to be determined by some 
underlying microscopic theory. In the case of the high-Tc cuprates, in partic- 
ular, there is no consensus about what would be such a theory, however it is 
generally accepted that superconductivity is constrained to the CuO'z planes. 

The complete lagrangian we will consider is given by 

^ = iV'aa ^i^aa + g (V'l\a ^^lla + V'lta V'ua) (V'2i6 V'lTb + V'Ub i^2^b) , (l) 

where (? > is a constant that may depend on some external control param- 
eter, such as the pressure or the concentration of some dopant. Later on, we 
will make g = jr. 



3 



We use the following convention for the Dirac matrices: 

/ = 7V = ^^ 7V = o-^ (2) 

Observe that the interaction lagrangian contains four terms, describing the 
various possible BCS interactions in which a Cooper pair would form between 
electrons with opposite spins, belonging to different Fermi points but in the 
same plane. Nevertheless, the interaction of electrons belonging to different 
planes is allowed. 

Furthermore, we can motivate our model out of a microscopic description of 
the high-Tc cuprates by noting that t-J model calculations of the electron 
spectral function indicate the emergence of small pockets of Dirac electrons 
at low doping [17]. 



In addition to a complex valued 0(N) 
possesses the U(l) symmetry 

and the chiral U(l) symmetry 



symmetry [18], the lagrangian above 
- e-'VL ^ = 1,2 (3) 

'lp2aa e-'^1p2aa- (4) 



We shall see below that the former symmetry is spontaneously broken at 
T = 0. A model presenting the spontaneous breakdown of the latter has been 
studied in [19]. 

We now introduce a Hubbard-Stratonovitch complex scalar field a through 



2a^J P(7*P(7exp|-i J Sx-^ (jvj 

= y P(7*P(7exp |-iy d^X-^ [a* - g (V'J^^ V2ia + ii^a ^lia) 

In terms of this, we may write the partition function 



as 



where 



(5) 

(6) 
(7) 



4 



^From this we obtain the field equation for the auxihary field a: 



(8) 



(9) 



As we shall see, the vacuum expectation value of a will be taken as the order 
parameter for the superconducting phase. 

We will now integrate over the fermionic fields. In order to do that, we intro- 
duce the Nambu fermion field = {'^iia '^lia V'lj.o V'lj.a) ■ terms of this we 
can rewrite (8) as 



where the matrix A is given, in momentum space, by 



(10) 



A 



-ko k- -a\ 
k+ —ko —o 
-a* -ko -k. 

-a* -A;_ —ko j 



with k± — vy {k2 ± iki)- 

Integrating on the fermion fields and redefining the coupling constant as = 
4, we obtain 



where 



N 



err -iA^lnDet 



A\a\ 
A[Q\ 



(12) 



(13) 



The determinant of the matrix A is dety4.[cr] = (/cq — '^fI^P) 
the above expression becomes 



— cr 



, hence 



S. 



eff 



d^xi-^\af ] -i27VTrln 



1 + 



□ 



(14) 
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3 The Superconducting Transition at T = and Quantum Criti- 
cality 



Let us consider in this section the T — case. We shall see that a quantum 
phase transition occurs, connecting the superconducting and normal phases. 



3.1 The Renormalized Effective Potential 



Since we are considering the zero temperature case, the functional integral 
in (12) must be dominated by constant configurations of a, which minimize 
the effective potential per plane, Ve«, corresponding to (14). This is more 
conveniently evaluated in the euclidean space and is given by 



Kff(kl) 



A 



-2 




1 + 



+ VpA;^ 



(15) 



where, henceforth, by a we actually mean (0|(j|0). The above expression cor- 
responds to a mean field approximation. Conversely, this would be the leading 
order approximation in an 1/N expansion and would be the exact result for 
^ oo. 



The explicit form of the effective potential may be evaluated in this frame- 
work, by introducing a large momentum cutoff A/vf- The resulting expression, 
obtained from (15) is 



Kfr(kl) 



\a\ 



A \a\' 



+ 



(16) 



dada* 



Quartic fermionic theories in 2-l-lD have been shown to be renormalizable in 
& 1/N expansion [20]. In order to eliminate the divergent constant A from 
(16) we renormalize the coupling constant A using the usual renormalization 
condition [21] 

= ^ (17) 

|(T|=fT0 

where ctq is an arbitrary finite scale parameter, the renormalization point and 
Xr is the (finite), renormalized coupling constant. 

Inserting (16) in (17), we obtain 
where a = 27rvp. 
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Substituting this result in (16), we get the renormahzed effective potential per 
plane 



3.2 The Gap Equation and the Quantum Critical Point 



Let us study now the minima of the renormahzed effective potential per plane, 
Eq. (19). The first and second derivatives of K;ff,R with respect to |cr| are given, 
respectively, by 

v;i,K(kl) = 2H(^-^.M) (.0) 

and 

^-(H)^2(^-|f.^) (21) 

The ground state is determined by the solutions of V^g ^ = 0. This admits two 
solutions, which we call A. Notice that 

A=|(0|a|0)| (22) 

and (0|(7|0) is a complex order parameter for superconductivity. 

^Prom (20) we conclude that either A = or A 7^ 0, the nonzero solutions 
satisfying the gap equation 

l = ^f^-AV (23) 



a V 2 

Inserting the A = solution in (21), we get 

Ki,(A^0,^2(J--|^) (24) 

and we conclude that A = will be a minimum only for Xr < 2a /Sao = ^c- 
^From (23), conversely, we see that the A 7^ solution is given by 

Since Aq is positive semi-definite, the above expression will actually be a 
solution only for Ar > Ac. On the other hand, from (21), we immediately see 
that V^g (A) = 2Ao/a > for the solutions of the gap equation (23). 
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We can infer that the ground state of the system will be 







Xr < Ac 



(26) 



) 



Xr > Ac 



Expression (26) implies that the system undergoes a continuous quantum 
phase transition at the quantum critical point Ac = 47rf|/3cro. Since A is 
the modulus of the order parameter for superconductivity, we conclude that, 
for couplings below Ac the electronic system will be in the normal state, while 
for couplings above Ac, it will be in a superconducting one (for temperature ef- 
fects, see teh next section). The quantum critical point Ac, therefore, separates 
a normal from a superconducting phase. 



3.3 Renormalization Group Analysis 

Let us now apply renormalization group methods, in order to show that our 
results are completely independent of the arbitrary finite scale ctq introduced 
in our renormalization procedure. 

We start showing that the renormalized effective potential Ves,R, given by (19) 
satisfies a renormalization group equation. Indeed it is easy to see that 



with the /3-function given hy P — — 3(7o/2q;A|j = — A|j/Ac. This means that 
the renormalized effective potential does not depend on the renormalization 
point (To- A negative /5-function, on the other hand, implies that the theory 
is asymptotically free. Indeed, keeping Aq fixed and taking the limit ^ oo 
we see that A^j — > 0. 

One can also show that the gap Aq, given by (25) satisfies the same renormal- 
ization group equation being, therefore, also independent of cxq. Indeed, using 
(25) and the expression of the /5-function just found, we obtain 



Finally, solving the differential equation corresponding to the /3-function def- 




(27) 




(28) 
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inition, namely, 



= /3, (29) 



we obtain 

1 1 1 



Xnicr'o) Xcia'o) Xnia'^) XM 



, (30) 

o) 

where ~ ^ '^o arbitrary values of the renormaliza- 

tion point. The above equation, clearly shows that the superconducting gap 
Aq, given by (25) is independent of the scale ag. 

The theory predicts the existence of a quantum critical point Ac, separating 
two phases at T = 0: a normal (Aq = 0) and a superconducting one (Aq 7^ 
0). Nevertheless, the theory docs not predict the value of Ac- This has to 
be determined experimentally. The renormalization group analysis, however, 
guarantees that the physics of the quantum phase transition will not depend 
on the renormalization point ctq. 



3.4 Gaussian Quantum Fluctuations 



Let us now investigate whether the results we found in the saddle-point ap- 
proximation for the ground state of the fermionic system described by (1) are 
robust against quantum fluctuations at T = 0. Notice that for the actual ex- 
istence of the quantum phase transition it is necessary that the normal phase 
found in the saddle point approximation should not be removed by higher 
order corrections, as it happens, for instance, in the Coleman- Weinberg mech- 
anism [21]. For this it is crucial that the corrected gap equation admits a 
Aq = solution. 

Expanding the effective action (14) in (12) about a stationary point cr and 
retaining the gaussian quantum fluctuations rj, we obtain, after integrating 
over T] and rj*, 

4ff[(j]=5efr[(j]-TrlnM[(7], (31) 
where Ses [c] is given by (14) and 



M 



( s^s s^s \ 
SaSa* Sa"^ 
S^S 
(5ct*2 Sa*SaJ 



(32) 



The effective potential corresponding to this is given by 

v;ff(ki) = Kff(ki) + v(|<7|), 

where V^fr {\a\) is given by (16) and 



(33) 
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{27ryJ 27r 



A 



— a(a;, k) 



+ ln 



A 



a{uj,k) +2\a\%{u,k) 



}• 



with 



d^q f de 



a{u;,k,\a\) = J-^J 



(27ry J 27r[^;p|^ +i^] + k\ - i {6 + u) 



+ k 



(34) 



(35) 



and b{u!, k, \a\ 



d a 



We would like to stress at this point that, for the expansion about the saddle- 
point to make sense, we must have V (a) <S Ves (a) . This condition is actually 
usually satisfied because the expansion is a power series in h. Keeping this 
fact in mind, let us look for the minima of Ves (cr). From (34), we can infer 
that the first derivative of V (cr) is of the form 



V'(M) = M fi\a\), 



(36) 



where /(O) is a finite constant. Combining (36) with (20), we immediately 
conclude that Aq = is a solution of the corrected gap equation V^g (a) — 0. 
The second derivative of the corrected effective potential V^fr (cr), evaluated at 
this point is 



v;'4(A = o) = 2 



i__^)+V''(A..O). 



(37) 



where wc have used (24). Since yrr cc h?,it follows that the second term in the 

cff 

rhs of the above equation must be much smaller than the first one. Therefore 
the sign of t^s (c — 0) must be the same as the one of T^ff (cr = 0). We conclude 
that Ao = is indeed a minimum, even considering the quantum fiuctuations. 
For the same reason, the mean field solution Aq 7^ should not be cancelled 
by quantum corrections. Thus, we conclude that the phase structure we found 
at T = is robust against quantum fiuctuations. 



4 The Superconducting Transition at T 7^ 



4.1 The Gap Equation atTy^O 



We consider in this section the finite temperature effects in the superconduct- 
ing transition and in the order parameter A. The nonzero solutions for A 
at a finite temperature are supposed to hold a priori only in the — > 00 
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limit, because otherwise they are ruled out by the Coleman-Mcrmin-Wagner- 
Hohcnbcrg theorem [16]. This limit corresponds to a physical situation where 
the three-dimensionality of the system is explicitly taken in account. For finite 
values of and T 7^ 0, the situation is quite subtle. We discuss it in the next 
subsection. 

At T 7^ 0, the effective potential is no longer given by (15). It must be replaced 
by 



Kff(|a|,T) 



A 



2T 



E In 



1 + 



(38) 



where uin = {2n+l)TiT are the fermionic Matsubara frequencies corresponding 
to the functional integration over the electron fields. 

The finite temperature corrections do not alter the ultraviolet divergence 
structure of the theory, hence we may eliminate the divergences in the T ^ 
case through the same renormalization of the coupling constant A as in the 
zero temperature case, given by (18). 

In order to derive the gap equation, we consider the following condition, which 
must be satisfied by the order parameter at a finite temperature: 



K'ff(|<T|,r)==2|<7| 



A 



1 

2a Jo 



dx 



X + cr 



tanh 



Ix + a 
2T 



0. (39) 



This was obtained by taking the derivative of (38) with respect to \a\ and 
performing the Matsubara sum. In the above equation, A is the same cutoff 
used before. 



As in the T — case, this admits two solutions, either A(T') = or A(T') 7^ 
0. In the latter case, the superconducting order parameter satisfies the gap 
equation 

1 = - /^dytanh (— 
aJA ^ \2T, 



(40) 



Solving the integral and renormalizing the coupling A as in (18), we find a gen- 
eral expression for the superconducting gap as a function of the temperature, 
namely 



A(T) = 2rcosh-^ 



e 2T 



(41) 



where Aq is given by (26), see Fig. 1. From (41) we can verify that indeed 
A(T = 0) = Aq. Also from the above equation, we may determine the critical 
temperature Tc for which the superconducting gap vanishes. Using the fact 
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that A(Tc) = 0, we readily find from (41) 



= (42) 
' 21n2 ^ ' 

This relation has been previously found in systems with a natural cutoff in the 
weak coupling regime [12] (see next section). It has also appeared in systems 
with dynamical mass generation in 2+ ID [22,23,24,25]. 

In Fig. 2, using (26) and (42), we display Tc as a function of the coupling 
constant. This qualitatively reproduces the superconducting phase transition 
of the high-Tc cuprates in the underdoped region. Since our theory describes 

the generic superconducting interaction of two-dimensional Dirac electrons, 
we may see this result as an indication of the possible relevance of this type 
of electrons in the high-Tc mechanism. 

In terms of the critical temperature, we may also express the gap as 

A(r) = 2rcosh-^ 2(^-^) . (43) 

Near this yields 

A(r) '"^'^ 2 V2h^ (l - ^) ' , (44) 

which presents the typical mean field critical exponent 1/2. 

We would like to remark, finally, that both the gap A(T) and the critical 
temperature do not depend on the arbitrary renormalization point ctq. 



4..2 Dynamical Gap Generation versus Spontaneous Symmetry Breaking 

The well-known Coleman-Mermin-Wagner-Hohenberg theorem [16] forbids the 

occurrence of spontaneous breakdown of a continuous symmetry at a nonzero 
temperature for systems in two spatial dimensions. Our superconducting order 
parameter is complex and given by 

(0|(7|0) = Ae^^ , (45) 

where A is the gap. Given the form of the field a, (9), we infer from (45) 
that a nonzero value for the gap would imply, in principle, the spontaneous 
breakdown of the U(l) symmetry, (3). For N oo, the system is effectively 
three-dimensional and the occurrence of a nonzero gap A(T) as determined in 
this section is is not in conflict with the theorem. For T = 0, either for flnite 
or infinite N, also we have a non- vanishing superconducting gap, which was 
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studied in section 3. In both cases the Coleman-Mermin-Wagner-Hohenberg 
theorem does not apply and a nonzero gap leads to a nonzero order parameter 
according to (45). 

In the case of finite N and T 7^ 0, in order to comply with the theorem and 
yet having a superconducting phase, we may invoke the mechanism proposed 
by Witten [26], by means of which we may have dynamical generation of a 
superconducting gap without the corresponding U(l) symmetry breakdown. 
It goes as follows. Whenever the gap is nonzero, according to (45), we must 
shift the field a as 

a^a- Ae'^ (46) 

in (8). This will produce an extra term in the effective lagrangian (8). In terms 
of new fermion fields, defined as 



'ipiaa = e '2 Vicra, (47) 



the extra term in (8) reads 



A 



(V^lTa i^lla + ^2Ta ^Ua) + (V^2ib i^l^b + V^Ub V^2Tb)] • (48) 



This is an explicit superconducting gap term that will make 

{0\i^ha4la + 4^ai^lla\0)^0. (49) 

Since the U(l) symmetry acts as ipi„a e^'^'ipiaa and 6 ^ 6 -\- 2uj we imme- 
diately see that the field ipi^^a is invariant under U(l) rotations and therefore 
the non-vanishing expectation value above does not imply spontaneous break- 
down of the U(l) symmetry (notice that the chiral U(l) symmetry (4) is also 
unbroken). Thus, we can have dynamical generation of a superconducting gap 
without the associated spontaneous symmetry breaking [26]. 

We must examine the thermodynamic conditions for the occurrence of this 
situation. This has been done in detail for the case of the Gross-Neveu model in 
2+ ID [25] and also in the case of the semimctal-cxcitonic insulator transition 
that occurs in layered materials [29], both related to the potential spontaneous 
breakdown of the chiral symmetry. The results of these analysis also apply 
here. 



The basic point is that, in order to check whether the order parameter (45) 
is zero or not, we must analyze the thermodynamics of the phase 6 of the 
superconducting order parameter. It turns out that this phase decouples and 
suffers a Kosterlitz-Thouless [27] transition at a temperature T^t- For temper- 
atures above Tkt there is no phase coherence and the superconducting order 
parameter vanishes because then {cos 6) = (sinO) = (even though A may 
be different from zero). Below Tkt there is a phase ordering and there will 
be a nonvanishing gap provided the condition T < T(. is also met (otherwise 
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A = 0). As it is, Tkt < Tc [25] and, therefore, the actual super conduct- 



clearly indicates that in spite of the fact that we may have a superconducting 
gap at a finite temperature in two-dimensional space, only in a really three- 
dimensional system we will have phase coherence developing at the same time 
that the modulus of the order parameter becomes nonzero, as determined by 
the gap equations. 

It has been speculated [28] that the above scenario could provide a framework 
for explaining the pseudogap transition that precedes the superconducting 
transition in high-Tc cuprates in the underdoped region. Our model provides 
a concrete realization of this mechanism. 



5 Systems with a Physical Cutoff A 

5.1 Physical Cutoff 

When considering applications in condensed matter systems, one usually finds 
a natural energy cutoff A (momentum cutoff A/vp). The Debye frequency 
(energy) is an example, in the case of conventional BCS superconductors. In 
this case, no renormalization is needed and the coupling constant A is the 
physical one. We investigate in this section the modifications that will occur 
in the superconducting electronic system under consideration when there is a 
physical cutoff in energy or momentum. 

The two-body interaction in this case, instead of being a delta function leading 
to the local interaction in (1), is given, in terms of the momentum cutoff, by 



ing transition occurs at T^t- It can be shown that Tkt 



N^oo 



Te [25]. This 



- A 



\k\ < A/vf 



V{k) = < 



(50) 







1^1 > A/vp 



In coordinate space, this corresponds to the interaction potential 




(51) 



where Ji is a Bessel function. 
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5.2 TheT^O Case 



We must now evaluate (15) with a finite physical momentum cutoff A/vp. This 
yields 



Kfr(kl) = 



(kr+A')' 



(52) 



A 3a 

We would like to stress that we are not assuming that A is large compared to 
I (J I, rather, we are considering a completely arbitrary finite cutoff A. This is not 
the situation usually found in conventional BCS superconductors. However, 
it is likely to be found in nonconventional ones such as high-Tc cuprates. For 
A ^ \a\ (52) would reproduce (16). 



The solutions of 



2\a\ 



|a|2 + A2) 



2^l 



a 



+ — 
a 



(53) 



will give the gap in the present case. This admits two solutions, namely, Aq = 
or 

A -^\—-—' 



(54) 



The second derivative of the potential, evaluated at Aq = is 

A 



yi(A. = o) = i-- 



(88) 



and we conclude that Aq = is a solution only for A < Ac, with Ac = a/A. 
Conversely, the second derivative of (52) evaluated at Aq 7^ ( given by (54)) 
is positive for A > Ac. As a consequence the gap now will be given by 



ctA 
2 







A? 



A < Ac 



(56) 



A2 



A > Ac 



which should be compared with (26), see Fig. 3. Again we find a quantum 

phase transition, now at the critical coupling Ac = a/A. Observe now that, 
contrary to the local case, since A is a physical parameter, the value of the 
quantum critical point Ac is predicted by the theory. 

Observe that, for A > Ac 



Ao 



a 



,Ac 



1^ 



'A + Ac^ 
, 2Ac , 



(57) 
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and in the region where A > Ac, we have Aq — Aq, with Aq given by (26). 
This coincides with the regime where A ^ Aq and occurs near the quantum 
critical point. Also, from (57) we see that the quantum phase transition is in 
the same universality class as the one found in Sect. 2, as it should. 



5.3 TheT^O Case 



Let us consider now the case of systems with a physical cutoff at a finite 
temperature. In this case the gap equation (58) must be replaced by 



(58) 



Solving the integral, we find an implicit equation for the superconducting gap 
in the presence of a physical cutoff A, at an arbitrary temperature, namely 



A(T) = 2Tcosh-^ 



e 2TA cosh 



A2(r) + A^ 

27^ 



2 \ 2 



(59) 



Using the fact that Aq satisfies (53), we may verify that indeed 



A(T) 2rcosh-^ 



e 2T 



(60) 



where Aq is given by (56). 



Let us now determine the critical temperature Tc, for the onset of supercon- 
ductivity in the presence of a physical cutoff. We must have A(Tc) = and 
therefore, from (59) 

cosh (^— ) = . (61) 



This equation yields the following relation between Tc and the zero tempera- 
ture gap, (56) 

Ao= f^l2Tcln[^-^ , (62) 



2A. 



l + e 



where Xc — a/A. 



A particular regime that is frequently studied is the one where the cutoff is 
large compared to the critical temperature and to Aq, namely, when A ^ 
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Tc and A > Ac- Observe that this last condition guarantees that A ^ Aq, 
according to (56) and (57). Since the previous relations imply 

7; « A = ^ ~ -, (63) 

Ac A 

we may infer that the conditions above hold in the weak coupling regime. 
In this case, (62) becomes simply Aq = 21n2 [12]. Also, according to 
(56) and (57) Aq — > Aq in this regime. Thus wc recover (42), the relation 
previously found between Tg and the zero temperature gap in systems without 
a natural cutoff . This relation gives the ratio Aq/Tc ~ 1.39, which should be 
compared with the corresponding ratio in the BCS theory for conventional 
superconductors, namely, 1.76, which is also derived in the limit A/Tc ^ 1. 

In the weak coupling regime, given by (63), we also recover expressions (43) 
and (44), for the superconducting gap as a function of the temperature. The 
pre-factor in the latter expression, describing the behavior of the gap around 
Tc is 2.36, whereas the corresponding value in BCS theory is 3.06. 

We solve (59) numerically for A, using different values for the ratio A/ Ac and 
display the result in Fig. 4. In this, we may observe that indeed in the weak 
coupling regime the ratio Aq/Tc approaches the result given by (42). On the 
other hand, as the coupling increases, we see that this ratio surpasses the 
value obtained in BCS theory and approaches the values obtained in strongly 
coupled systems. 

There are condensed matter systems for which the weak coupling condition 
(63) is not valid. Should one of such systems contain Dirac fermions, we should 
use (56) and (59) for the superconducting gap, respectively at T = and 
T 7^ 0. The critical temperature, by its turn, would be given by (62). 

It is remarkable that the expression we find for the superconducting gap in the 
case of Dirac fermions, strongly differs from the one obtained in BCS theory in 
any dimensions. There the gap has an exponential dependence on the inverse 
of the coupling A. Here, in spite of still being non-analytical in the coupling 
constant A, the gap has a power-law dependence on it and vanishes below a 
critical value at zero temperature. The different behavior can be traced back 
to the fact that in the case of Dirac fermions the density of states vanishes at 
the Fermi points. In the BCS case, however, we have a Fermi surface with a 
non- vanishing density of states at the Fermi level and the momentum integrals 
used for obtaining the effective potential may be evaluated as 

p^^^(*)X/e (64) 

where N(eF) is the density of states at the Fermi level. This leads to a gap 

1 

proportional to e ^^^f^-^ . 
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6 Conclusions 



Our results highlight the qualitative and quantitative differences existing be- 
tween Dirac electrons - for which the Fermi surface reduces to a point - and 
quasi- free electrons, having a dispersion relation e{k) — \k\'^/2m*, which leads 
to a Fermi surface formation. The properties of the former are analyzed when 
the conditions are such that a superconducting interaction is present in a 
quasi-two-dimensional system consisting of a stack of N planes. One of the 
most striking differences between the two electron types, is the polynomial, 
rather than exponential gap dependence on the inverse coupling constant, in 
the case of Dirac electrons. This leads to a quantum phase transition separat- 
ing a normal from a superconducting phase for a critical value of the coupling, 
at T = 0. This transition is rather similar to the one occurring in the Nonlin- 
ear Sigma model in 2-l-lD, where the spin stiffness has an expression identical 
to (26) [30]. Because of this quantum phase transition, our model reproduces 
qualitatively the superconducting transition in high-Tc cuprates in the under- 
doped regime as one can infer from Fig. 2. This seems to indicate that Dirac 
electrons may have an important role in the mechanism of high-Tc supercon- 
ductivity. 

The ratio between the zero temperature gap to T^, as well as the behavior of 
the gap around Tg, are results also significantly different for the case of Dirac 
electrons. It would be very interesting to compare our results with experi- 
mental measurements of these quantities in quasi-two-dimensional materials 
containing Dirac electrons. 

An interesting outcome of this work is the analysis the phase structure gen- 
erated by the model in the presence of a natural physical cutoff, specially the 
regime where the cutoff is of the same order of the gap and we are away from 
the quantum critical point. In this case the (naturally) cutoff theory yields 
results that are strongly different from the ones derived in the weak coupling 
regime where the cutoff is much larger than the gap and the system is close 
to the quantum critical point. The strong coupling regime, in particular, is 
probably relevant for the high-Tc cuprates. 



We would like to thank S.Sachdev and A.H.Castro Neto, for very helpful 
comments and conversations. 

This work has been supported in part by CNPq and FAPERJ. ECM has been 
partially supported by CNPq. LHCMN has been supported by CNPq. 



18 



References 



[I] S.H.Simon and P.A.Lee, Phys. Rev. Lett. 78 (1997) 1548; A.C.Durst and 
P.A.Lee, Phys. Rev. B 62 (2000) 1270 

[2] E.J.Ferrer, V.P.Gusynin and V. de la Incera, Mod.Phys.Lett. B 16 (2002) 107 

[3] I.F.Herbut, Phys. Rev. Lett. 88 (2002) 047006 

[4] M.Franz and Z.Tesanovic, Phys. Rev. Lett. 84 (2000) 554, Phys. Rev. Lett. 87 
(2001) 257003 

[5] P.W.Anderson, l cond-mat/9812063| (unpubhshed) 

[6] G.Semenoff, Phys. Rev. Lett. 53 (1984) 2449 

[7] F.D.M.Haldane, Phys. Rev. Lett. 61 (1988) 2015 

[8] J.Gonzalez, F.Guinea and M.A.H.Vozmediano, Phys. Rev. Lett. 69 (1992) 172; 
Nucl. Phys. B406 (1993) 771; Nucl. Phys. B424 (1994) 595; Phys. Rev. Lett. 
77 (1996) 3589; Phys. Rev. B 59 (1999) 2474; Phys. Rev. B 63 (2001) 134421 

[9] N.M.R.Peres, F.Guinea and A.H.Castro Neto, |cond-mat/0506709l 

[10] Y.Zhang et al., Nature 438 (2005) 201 

[II] K.S.Novoselov et al.. Nature 438 (2005) 197 

[12] A.H.Castro Neto, Phys. Rev. Lett. 86 (2001) 4382 

[13] B.Uchoa, A.H.Castro Neto and G.G.Cabrera, Phys. Rev. B 69 (2004) 144512; 
B.Uchoa, G.G.Cabrera and A.H.Castro Neto, Phys. Rev. B 71 (2005) 184509 

[14] L.Balents, M.P.A.Fisher and C.Nayak, Int. J. Mod. Phys. B 12 (1998) 1033 

[15] L.Balents and M.P.A.Fisher, Phys. Rev. B 55 (1997) R 11973 

[16] N.D.Mermin and H.Wagner, Phys. Rev. Lett. 17 (1966) 1133; P.C.Hohenberg, 
Phys. Rev. 158 (1967) 383; S.Coleman, Commun. Math. Phys. 31 (1973) 259 

[17] X-G.Wen and P.A.Lee, Phys. Rev. Lett. 76 (1996) 503 

[18] E.C.Marino and M.J.Martins, Phys. Rev. D 33 (1996) 3121 

[19] G.W.Semenoff and L.C.R.Wijewardhana, Phys. Rev. Lett. 63 (1989) 2633 

[20] B.Rosenstein, B.J.Warr and S.H.Park, Phys. Rev. Lett. 62 (1989) 1433 

[21] S.Coleman and E.Weinberg, Phys. Rev. D 7 (1973) 1888 

[22] A.Okopinska, Phys. Rev. D38 (1988) 2507 

[23] B.Rosenstein, B.J.Warr and S.H.Park, Phys. Rev. D39 (1989) 3088 

[24] T.Appelquist and M.Schwetz, Phys. Lett. B491 (2000) 367 



19 



[25] E.Babaev, Phys. Lett. B497 (2001) 323 

[26] E.Witten, Nucl. Phys. B145 (1978) 110 

[27] V.L.Berezinskii, Zh. Eksp. Teor. Fiz. 59 (1970) 907; J.Kosterlitz and 
D.Thouless, J. Phys. C6 (1973) 1181 

[28] V.M.Loktev, R.M.Quick and S.G.Sharapov, Phys. Rep. 349 (2001) 1 

[29] D.V.Khveshchenko and H.Leal, Nucl. Phys. B687 [FS] (2004) 323 

[30] S.Chakravarty, B.I.Halperin and D.Nelson, Phys. Rev. B 39 (1989) 2344 



20 




22 




23 




— vx = o.i 

1JX = 0.5 

- I A = 0.99 



0.0 



T/T. 



Fig. 4. The superconducting gap A(T) divided by Tc as a function of the normahzed 
temperature T/Tc for several values of Ac/A. The two arrows indicate the values for 
Aq/Tc given by the BCS theory and the A ~ Ac case (Ac = a/ A). 
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